AMENDMENT UNDER 37 C.F.R. §1.111 
U.S. Appln. No. 10/766,011 
Attorney Docket: Q79665 

AMENDMENTS TO THE CLAIMS 

This listing of claims will replace all prior versions and listings of claims in the 
application: 

LISTING OF CLAIMS: 

1 . (currently amended): A gene expression state estimating system for estimating the 

probability of gene expression in each channel, the system including an input device for sending 
gene expression level data, a program-controlled data analyzer, and an output device, wherein 
said data analyzer comprises 

distributed parameter estimating means for estimating distributed parameters of a mixed 
normal distribution shown in the following equation (25) using the gene expression level data 
from said input device, and sending the estimated distributed parameters: 

(i-^("-ftko 2 )+^(«-//,k, 2 ) 

(1 - l) q> (u - up | Uq) + gqi (u - | u^) ^ 

where </> (u - | <r 0 2 ) -${tt — ju Q | cr 0 2 ) represents the density function of a one- 
dimensional normal distribution with average j±o and variance <r 0 2 , (/> (u - ju x \ a] ) represents the 
density function of a one-dimensional normal distribution with average \i\ and variance a] a 
(M\ > °l ) (/*<)> °"o ) an< 3 (M\ > °\ ) are average and variance parameters of first and second 

components, respectively, u represents data points corresponding to sum of gene expression level 
from two channels and £ is the mixing ratio, with the assumption that 

Ho < Ha. °l > 0, o\ > 0, 0 < 5 < l issatisfied> 
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mixing ratio parameter estimating means for estimating a mixing ratio parameter of the 
mixed normal distribution using the gene expression level data sent from said input device and 
the distributed parameters sent from said distributed parameter estimating means, and sending 
the estimated mixing ratio parameter, and 

posterior probability calculating means for calculating the posterior probability of the 
expression state of each gene in each channel using the gene expression level data, the estimated 
distributed parameters and mixing ratio parameter, and sending the calculated posterior 
probability to said output device. 

2. (original): The system according to claim 1 wherein said distributed parameter 

estimating means estimates the mixing ratio (£), average (jxo, |Ui), and variance v °' i; by 
applying the mixed normal distribution of two components to data on the sum of the amounts of 
expression of genes located in a region where the difference of gene expression levels X and Y 
of two channels is near 0. 

3. (previously presented): The system according to claim 2 wherein when the median 
value of the absolute difference |vi|(i = l,...,n) of the gene expression levels X and Y is c M , the 

data on the amounts of gene expression is shown by {u } \ | Vj | < c M , i=l n} , and v represents 

data points corresponding to difference of the gene expression level of the two channels. 
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4. (original): The system according to claim 3 wherein said distributed parameter 

c *\ *\ 2 ^2 

estimating means performs estimation by the use of the estimated a °' Pl# 01 according 

2 2 

to the following equations (26), (27), (28), and (29) to estimate p ' ° e ' CT P' X ; 

a = (Pi-p 0 )/2 (26) 

-2 _ 1 v 1 -.2 

(27) 



o" = — n — n- y Vi 

2 KI is % 



A 2 1 £2 1 -2 

4 2 (28) 



X = (log 



Oi — CJ-> 
1 + ? 



. ^ 2 
4]i 



(29) 

where N 0 denotes an index set of data values that satisfies 1 e ^ I u i < £o) anc j 
II N °H denotes the number of elements. 

5. (previously presented): The system according to claim 4 wherein said mixing ratio 
parameter estimating means estimates the mixing ratio parameter p=(poo> Pio, Poi, Pn) (where p 0 o 
denotes a mixing ratio when no gene is being expressed in both cells 1 and 2, pi i denotes a 
mixing ratio when any gene is being expressed in both cells 1 and 2, p !0 denotes a mixing ratio 
when a gene is being expressed in cell 1 but not in cell 2, and p 0 i represents a mixing ratio when 

q * ^ ^2 a 2 v 

a gene is being expressed in cell 2 but not in cell 1) using " ^' ' ° c ' °P ' given from said 
distributed parameter estimating means by applying a two-variable mixed normal distribution of 
four components shown in the following equation (30) to gene expression level data {(u j? Vj)|i=l, 
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n} sent from said input device wherein v represents data points corresponding to difference 

of the gene expression level of the two channels. 

Poo9oo< u ' v | S) + Pi 0 gio(u, v | §) + Poig 0 i(u, v | §) + p n g n (u, v | 9) 

- p 00 q>(u | 4o£ + 2o|)<p(v | 2a 2 e ) + p 10 <p 2 (u - p, v - p | E 10 ) 

+ Poi<P2< u - p, v + p | E 01 + p lx <p(u - 2p | 4p 2 (e X - 1) 
+ 4a 2 , + 2o 2 ) 9 (v | 2a 2 ) 

where the above equation satisfies the relationships shown in the following equations 
(31) and (32). 



=10 = 



'$ 2 (e* 2 - 1) + 45^ + 2o 2 YL 2 (e p - 1) 



£ 2 (e x 



E oi - 



p 2 (e x - 1) + 43^ + 2o 2 - fi 2 (e x ' - 1) 
- p 2 (e x2 - 1) fi 2 (e x2 - 1) + 2o 2 



1) 



,-2 



p 2 (e x - 1) + 2o\j 



(31) 



J 



(32) 



6. (previously presented): The system according to claim 5 wherein said posterior 
probability calculating means calculates the posterior probability of expression of any gene in 
cell 1 and cell 2 for each pair (u, v) of the gene expression level data sent from said input device 

according to the following equation (33) (where f (u# v I P' ®) i s given by the following 
equation (34)) and the following equation (35) (where gene expression level data are shown in 
the equations (X = Tia + p + 6|,Y = x 2 a + p + e 2 ), and tj, t 2 take either 1 or 0, which represents 
the presence or absence of true gene expression in each cell). 

Pr (t x = 1 | p, §) = Pio g io (Uy v 1 g > + PiiSTii^ v I 6> 

f(u, v | p, §) m >. 
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f (u, v | p, e) 



(j,k)e{0,l} 2 



(34) 



Pr (t 2 = 1 | p, 8) = 



Poigpi^ v I §) + Pngn(vw v | 9) 
f(u, v l p, e> 



(35) 



7. (previously presented): The system according to claim 5 wherein said posterior 
probability calculating means calculates the posterior probability indicating an event of 
differential expression between cell 1 and cell 2 according to the following equation (36) 



where gene expression level data are shown in the equations (X = Tiot + |3 + ei,Y = 12a + 
P + £2), and t|, i2 take either 1 or 0, which represents the presence or absence of true gene 
expression in each cell. 

8. (original): The system according to claim 6 or 7 wherein said posterior probability 
calculating means judges whether calculations of posterior probabilities of gene expression have 
been made for all the pairs (u, v) of the gene expression level data, and when all the calculations 
have been completed, it ends the process, while when all the calculations have not been 
completed yet, it calculates the posterior probability related to the next gene, such that the 
calculated posterior probabilities of gene expression in each channel are sent to said output 
device, and 

said output device displays the posterior probabilities of gene expression in each channel. 



Pr (i! * t 2 I p, 9) = 



PipgiO^ v I 9) + p 01 g 01 (u, v I 6) 
f (u, v I p, 0) 



(36) 
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9. (original): A gene expression state estimating method for estimating the probability of 
gene expression in each channel based on gene expression level data, comprising the steps of: 

estimating distributed parameters of a mixed normal distribution shown in the following 
equation (37) using the gene expression level data and sending the estimated distributed 
parameters: 

(1 - 5) cp (u - p 0 | al) + 5<p (u - Pi | ol) (3?) 

where I ° ) represents the density function of a one-dimensional normal 

distribution with average 0 and variance a 2 , and are average and variance 

parameters of first and second components, respectively, and ^ is the mixing ratio, with the 

assumption that *° < Pi, °o > 0, o* >0, 0 < * < 1 is satisfieds 

estimating a mixing ratio parameter of the mixed normal distribution using the gene 
expression level data and the estimated distributed parameters, and sending the estimated mixing 
ratio parameter, and 

calculating the posterior probability of the expression state of each gene in each channel 
using the gene expression level data, the estimated distributed parameters, and the estimated 
mixing ratio parameter, and sending the calculated posterior probability. 

10. (original): The method according to claim 9 wherein said step of estimating the 
distributed parameters further comprises a step of estimating the mixing ratio (£), average (jio, 

(a 2 a 2 ) 

^i), and variance 0 ' 1 by applying the mixed normal distribution of two components to data 
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on the sum of the amounts of expression of genes located in a region where the difference of 
gene expression levels X and Y of two channels is near 0. 

1 1 . (original): The method according to claim 10 wherein when the median value of the 
absolute difference |vj| (i = l v ..,n) of the gene expression levels X and Y is c M , the data on the 
sum of the amounts of expression of genes is shown by {uj 1 1 Vi | < cm, i=l,...,n}. 

12. (original): The method according to claim 1 1 wherein the estimation is performed in 
said step of estimating the distributed parameters by the use of the estimated Po ' ° 0 ' Pl ' 0l 

2 2 

according to the following equations (38), (39), (40), and (41) to estimate p ' ° 8 ' °P f X ; 



(38) 



1 




i€N 0 



(39) 





(40) 




(41) 



where N 0 denotes an index set of data values that satisfies 1 e ^ I u i < £o) anc j 



denotes the number of elements. 
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13. (original): The method according to claim 12 wherein the estimation is performed in 
said step of estimating the mixing ratio parameter p=(poo, Pio, Poi, Pn) (where poo denotes a 
mixing ratio when no gene is being expressed in both cells 1 and 2, pi i denotes a mixing ratio 
when any gene is being expressed in both cells 1 and 2, pi 0 denotes a mixing ratio when a gene is 
being expressed in cell 1 but not in cell 2, and g 0 i represents a mixing ratio when a gene is being 

0 — (u A \ 

expressed in cell 2 but not in cell 1) using H/ 9 e ' P sent from said step of 

estimating the distributed parameters by applying a two-variable mixed normal distribution of 
four components shown in the following equation (42) to the sent gene expression level data {(m, 
Vi)|i=l v ..,n} 

Poo9oo< u ' v I §) + Pi 0 gio< u ' v | 8) + Poigoi(u, v | §) + p 11 g 11 (u, v | e) 

= Poo9(u | 4ag + 2o*)<p(v | 2b\) + p 10 <P 2 <u - p, v - p | E 10 ) 

+ Poi92< u - P, v + p | E 01 + p xl 9(u - 2£ | 4p 2 (e X - 1) 
+ iol + 2S2 )( p(v | 2o 2 ) 

where the above equation satisfies the relationships shown in the following equations 
(43) and (44). 



E 10 - 



p*(e" - 1) + 4a 2 , + 1o% p'(e A - 1) 



~2 



p 2 (e x ' - 1) 



p 2 (e x - 1) + 2o\ ) 



(43) 



'01 - 



V(e x2 - 1) + 4o|. + 2o 2 - p 2 (e^ - 1) 



~ 2 , A 



- p (e" - 1) 



p 2 (e X - 1) + 2o 2 ) 



(44) 



14. (original): The method according to claim 13 wherein the calculation of the posterior 
probability of expression of any gene in cell 1 and cell 2 is made in said step of calculating the 
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posterior probability for each pair (u, v) of the sent gene expression level data according to the 

following equation (45) (where f (u * v I P' e > \ s g i ve n by the following equation (46)) and the 
following equation (47) (where ti, x 2 take either 1 or 0, which represents the presence or absence 
of true gene expression in each cell). 

Pr (ti = 1 | p, §) = Pio g io (u * v I j? + Piifii^ v I §> 

f (u, v | p, §) (45 ^ 



f(u, v|p f §)= £ Pjk9jk(U/ v | 8) 

(j,k)e{0,l} 2 (46) 



Pr (x 2 = 1 | p, §) = Poi g oi (u# v 1 + Pngii (u * v I gj 

f (u, v | p, 6) ( 47 ) 



15. (original): The method according to claim 13 wherein the calculation of the posterior 
probability indicating an event of differential expression between cell 1 and cell 2 is made in said 
step of calculating the posterior probability according to the following equation (48) 

Pr(Ti * T2 |£,§) = PiQgiO^ v I S) + P 01 g 01 (u, v 1 §) 

f (u, v | p, 0) ^ 48 ^ 

where tj, t2 take either 1 or 0, which represents the presence or absence of true gene 
expression in each cell. 

16. (original): The method according to claim 14 or 15 wherein it is judged in said step 
of calculating the posterior probability whether calculations of posterior probabilities of gene 
expression have been made for all the pairs (u, v) of the gene expression level data, and when all 
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the calculations have been completed, the process is ended, while when all the calculations have 
not been completed yet, the posterior probability related to the next gene is calculated. 



17. (previously presented): A gene expression state estimating program for estimating 
the probability of gene expression in each channel based on gene expression level data, the 
program instructing a computer to execute the steps of: 

estimating distributed parameters of a mixed normal distribution shown in the following 
equation (49) using the gene expression level data and sending the estimated distributed 
parameters: 

(1 - 5) <p (u - p 0 | a?) + Sep (u - Pl | ol) (49) 

i 

where <f>(u-ju 0 \ cr 0 2 ) represents the density function of a one-dimensional normal 
distribution with average \x 0 and variance cr 2 , </> (u - ju } | cr 2 ) represents the density function of 

a one-dimensional normal distribution with average \i\ and variance a] , and 

(Pi/ <*i> are average and variance parameters of first and second components, respectively, u 
represents data points corresponding to sum of gene expression level from two channels and \ is 

2 2 

the mixing ratio, with the assumption that Po < Pl < a ° > 0/ 01 > 0/ 0 < 5 < 1 is 
satisfied, 

estimating a mixing ratio parameter of the mixed normal distribution using the gene 
expression level data and the estimated distributed parameters, and sending the estimated mixing 
ratio parameter, and 
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calculating the posterior probability of the expression state of each gene in each channel 

using the gene expression level data, the estimated distributed parameters, and the estimated 

mixing ratio parameter, and sending the calculated posterior probability. 

18. (previously presented): A computer-readable recording medium storing a gene 
expression state estimating program for estimating the probability of gene expression in each 
channel based on gene expression level data, the program instructing a computer to execute the 
steps of: 

estimating distributed parameters of a mixed normal distribution shown in the following 
equation (50) using the gene expression level data and sending the estimated distributed 
parameters: 

(1 - I) cp (u - p 0 | ag) + 59 (u - Hi | a?) (5Q) 
where (f> (u - ju 0 | cr 0 2 ) represents the density function of a one-dimensional normal 
distribution with average \x 0 and variance cr 0 2 , <f>{u - //, | a] ) represents the density function of a 

one-dimensional normal distribution with average \x\ and variance a] , and (//, ,cr 2 ) 

are average and variance parameters of first and second components, respectively, u represents 
data points corresponding to sum of gene expression level from two channels and £ is the mixing 

ratio, with the assumption that Po < Pi, «o > 0, oj > 0, 0 < * < 1 is satisfied? 

estimating a mixing ratio parameter of the mixed normal distribution using the gene 
expression level data and the estimated distributed parameters, and sending the estimated mixing 
ratio parameter, 
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calculating the posterior probability of the expression state of each gene in each channel 
using the gene expression level data, the estimated distributed parameters, and the estimated 
mixing ratio parameter, and sending the calculated posterior probability, and 

outputting the calculated posterior probability. 
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